Nonlocal Conductivity in Type-II Superconductors 

Chung- Yu Mou^, Rachel Wortis^'^, Alan T. Dorsey^ and David A. Huse^ 

^ Department of Physics, University of Virginia, 
McCormick Road, Charlottesville, Virginia 22901 
■ . ^ Department of Physics, University of Illinois at Urbana- Champaign, 

Q>^ . 1110 West Green Street, Urbana, Illinois 61801 

j^ . ^AT&T Bell Laboratories, Murray Hill, NJ 07974 

O' 
^ . (February 1, 2008) 

(N 

> 

'^ ■ Abstract 



0^ ■ Multiterminal transport measurements on YBCO crystals in the vortex 
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liquid regime have shown nonlocal conductivity on length scales up to 50 



1-^ ■ microns. Motivated by these results we explore the wavevector (k) depen- 

a: 

O ■ dence of the dc conductivity tensor, a^^Ck), in the Meissner, vortex lattice, 



and disordered phases of a type-II superconductor. Our results are based on 
time-dependent Ginzburg-Landau (TDGL) theory and on phenomenological 
arguments. We find four qualitatively different types of behavior. First, in 
the Meissner phase, the conductivity is infinite at k = and is a continuous 
function of k, monotonically decreasing with increasing k. Second, in the 
vortex lattice phase, in the absence of pinning, the conductivity is finite (due 
to fiux fiow) at /c = 0; it is discontinuous there and remains qualitatively 
like the Meissner phase for A; > 0. Third, in the vortex liquid regime in a 
magnetic field and at low temperature, the conductivity is finite, smooth and 
non-monotonic, first increasing with k at small k and then decreasing at larger 
k. This third behavior is expected to apply at temperatures just above the 

1 



melting transition of the vortex lattice, where the vortex liquid shows strong 
short-range order and a large viscosity. Finally, at higher temperatures in the 
disordered phase, the conductivity is finite, smooth and again monotonically 
decreasing with k. This last, monotonic behavior applies in zero magnetic 
field for the entire disordered phase, i.e. at all temperatures above Tc, while 
in a field the non-monotonic behavior may occur in a low temperature portion 
of the disordered phase. 
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I. INTRODUCTION 

In this paper we explore the nonlocal dc electrical transport properties of type-II super- 
conductors. What is meant by nonlocal here? The general expression connecting the local 
current density in a material, J, to the local electric field, E, in the linear (Ohmic) regime is 

J,{r) = Ja^,{r,r')E,{r')dr'. (1.1) 

If the conductivity o"(r, r') is nonvanishing for r 7^ r', then this is nonlocal; the current at 
r is determined by the field not only at r, but at all points r' where cr(r, r') 7^ 0. In a 
translationally invariant system, the nonlocal conductivity can only be a function of the 
difference (r — r'). Taking the Fourier transform one then obtains Jfj,(k) = crf^,y{'k)E,y{k). All 
materials exhibit nonlocal transport properties on some length scale. In normal metals, the 
nonlocal terms are significant only at length scales less than or of the order of the inelastic 
mean free path. In superconducting materials, however, as the transition is approached from 
above, the associated correlations can cause nonlocal effects to become important over much 
longer length scales. It is this phenomenon that we address in this paper. 

Vortices are important actors in the nonlocal electrical transport properties in type- 
II superconductors. This is sometimes more easily described in terms of the nonlocal dc 
resistivity:El 

E^{r) = Jp^,{r,r')Mr')dr'. (1.2) 

(Note that the nonlocal resistivity is a linear operator that acts on the full current pattern, 
J(r), and is, as usual, the inverse of the nonlocal conductivity operator.) A current J(r') 
pushes on the vortex segments at r' due to the Lorentz and Magnus forces. These vortex 
segments move, and, due to the continuity and entanglement of vortices as well as the 
repulsive and attractive forces between parallel and antiparallel (respectively) vortices, they 
cause vortex segments at r to move.0 This produces phase slip and electric fields at r. Thus 
vortex dynamics contribute significantly to the nonlocal resistivity, p^,^(r, r'). 



Recent multiterminal transport experiment slI on YBCO crystals in the vortex liquid 
regime obtain results that demonstrate nonlocal conductivity on length scales of at least 50 
microns. A phenomenological understanding of those experimental results was presentedlil 
based on a hydrodynamic description of the vortex liquid and its viscosity. In this paper we 
also examine the vortex lattice phase and extend the phenomenological study of the vortex 
liquid regime to low and zero magnetic field. We then calculate the nonlocal conductivity 
due to the Gaussian fluctuations above Tc within time-dependent Ginzburg-Landau (TDGL) 
theory. An outline of the paper is as follows: 

Section II considers an ideal, unpinned, defect-free vortex lattice. There the steady-state 
motion of the vortices is only a rigid-body (flux-flow) motion of the vortex lattice as a 
whole. The electric field due to this vortex motion is therefore determined only by the total 
force and torque on the entire vortex lattice. A dc current pattern with wavevector k ^ 
elastically distorts the lattice, but, in linear response, produces no steady motion of the 
vortices. Thus we find that in this ideal case the Ohmic flux-flow resistivity due to vortex 
motion is nonzero only at k = 0. In both the vortex lattice and Meissner phases there is an 
additional contribution to the resistivity that is proportional to k"^ at small k and does not 
arise from vortex motion.Q 

In Section III we discuss the phenomenology of the vortex liquid regime. Connectivity, 
entanglement and other interactions between nearly parallel vortices then give rise to the 
vortex-liquid viscosityo that reduces their motion in response to nonuniform dc currents, 
producing a resistivity that has a local maximum for uniform {k = 0) currents and is 
smaller for long-wavelength nonuniform {k ^ 0) currents, as in the case of the vortex lattice. 
However, at high enough temperature and low or zero magnetic field there will be thermally 
excited vortex lines and loops present with all orientations. Then the connectivity and 
attraction of nearly antiparallel vortices produces a nonlocal effect of the opposite sign, with 
the resistivity smallest at fc = 0. In the disordered phase, the conductivity and the resistivity 
may be expanded in powers of k; these arguments then indicate that the signs of the order 
k"^ terms will depend on temperature and magnetic field. 



Section IV sets up the TDGL calculation of the lowest-order (Gaussian approximation) 
fluctuation contribution to the nonlocal dc conductivity in the disordered phase. For a 
uniform {k = 0) current this Gaussian approximation gives the Aslamazov-Larkin fluctuation 
conductivity.!^ The calculation for zero magnetic fleld is carried out in Section V. Here the 
full fc-dependence can be obtained and the conductivity is found to be a monotonically 
decreasing function of k. Section VI obtains the fluctuation conductivity to order k'^ in a 
magnetic fleld. At this order, the sign of the fc- dependence remains the same as for zero 
magnetic fleld, at least for longitudinal electric flelds. We conclude with a brief summary 
and discussion in Section VII. Appendices A and B are devoted to some technical details of 
the TDGL calculation. 

II. VORTEX LATTICE 

Let us flrst consider the vortex lattice phase. In dc steady state the total time-averaged 
force on each portion of the vortex lattice, including drag forces, must vanish, since the 
lattice is either stationary or moving with a flnite steady-state velocity. In linear response 
to a dc current, the resulting steady-state local force balance equation for the vortex lattice 
contains four terms: one proportional to and perpendicular to the local current density J 
(arising from Lorentz and/or Magnus forces) ,0 one proportional to and perpendicular to the 
local vortex velocity v (arising from the Magnus force and/or a perpendicular drag force), 
one proportional to and parallel to the local vortex velocity (a standard drag force), and 
one arising from the local elastic distortion of the vortex lattice.111 Assuming the magnetic 
induction B = Bz, so that the vortices are parallel to the z-axis, and that the material is 
isotropic in the xy-p\a.ne, the steady state force balance equation is 

X B + 72nz X V — 7inv = 0, (2.1) 



where 



Helastic = 7^'^^ii-^){cL0^)kikj + 6ij[C(iQ{k)k]_ + C44(k)/c^]}Mj(k) (2.2) 



is the elastic energy of the vortex lattice.H u is the local displacement of the vortices away 
from an ideal, undistorted lattice, and c^, Cqq and C44 are the bulk, shear and tilt elastic 
moduli, respectively. Here v refers to the local velocity of the vortices in the steady state, 
averaged over time and over a length scale longer than the lattice spacing. We assume 
the lattice is dislocation free so that u is well defined, and that there are no vacancies 
or interstitials so that v is simply du/dt. The areal density n is related to the magnetic 
induction by n = B/(j)o, with 0o = h/2e the flux quantum. The drag coefficients 71 and 72 
are phenomenological parameters. 

When the current is spatially uniform, the resulting vortex velocities and displacements 
are also spatially uniform. The elastic term in ( p.l|) is then zero and we are left with a 
balance of the other three terms. Using the Josephson relation for the electric field produced 
by moving vortices,H 

E = -V X B, (2.3) 

one finds for the components of the conductivity tensor 

00-D 

- Jl. 

Gzz is infinite and the other components are zero. The effect is that when a spatially uniform 
current is applied, the vortex lattice as a whole is pushed across the sample, causing phase 
slip and therefore dissipation. This is the standard fiux-fiow resistivity. 

Now suppose the current is nonuniform. The elastic term is then nonzero because differ- 
ent parts of the vortex lattice are subject to different forces and this results in elastic strains. 
If the current produces a nonzero total force or torque on the vortex lattice, it will move 
as a (overdamped) rigid body in response, again exhibiting flux-flow resistivity. However, 
if the spatially averaged current (as well as the total torque applied to the vortex lattice 
by the current) vanishes, there is no steady-state motion of the vortex lattice. Instead the 
lattice is statically strained so that everywhere the elastic restoring force locally balances 



the force due to the nonuniform current. The vortices are not moving, so no electric field 
is being produced, although a current is flowing. Therefore, the resistivity due to vortex 
motion vanishes in linear response to a long-wavelength, static, nonuniform current pattern 
with zero spatial average. For example, if a dc current is applied in the x direction with 
variation in the x direction, then the vortices will be pushed along the y direction, producing 
static shear elastic distortions. Such an example is illustrated in Fig. 1. 

Here the dc resistivity due to vortex motion is completely nonlocal, since there is only 
rigid body motion of the vortex lattice as a whole, which is determined by the total force 
(and torque) applied to the entire vortex lattice by the current. Of course, this is a very 
idealized discussion, since we are neglecting pinning, defects in the vortex lattice (and the 
resulting plastic motion), and the boundaries of the sample. However, this shows that in this 
ideal case, the freezing of the vortices into a lattice is indeed a superconducting transition. 
Although the fiux-fiow resistivity to a uniform {k = 0) current does not vanish when the 
vortices freeze, the dc fiux-fiow resistivity to a long-wavelength nonuniform {k > 0) current 
does vanish. 

We have examined the response of the vortex lattice to a nonuniform electric field within 
the time-dependent Ginzburg-Landau (TDGL) equations in the absence of thermal noise, 
following, e.g., Troy and Dorsey.t3 The above phenomenology for the resistivity due to 
vortex motion is confirmed, but for nonuniform currents there is an additional contribution 
to the resistivity which does not arise from vortex motion. This dissipation in a static, but 
strained configuration arises directly from the gradient-squared term in the TDGL equations 
and gives a resistivity proportional to k"^ for small wavevector k in both the vortex lattice 
and Meissner phases. Q 

To summarize, the resistivity in the ideal vortex lattice phase is nonzero at A; = and 
proportional to k^ for small nonzero k. It is therefore discontinuous and non-monotonic 
in k a.t k = 0. We expect that this discontinuity and non-monotonicity is still present if 
the lattice contains a non-zero density of vacancies and interstitials: then the resistivity for 
k = is due to motion of the entire lattice, while for small nonzero k only the defects move. 



resulting in a lower resistivity. 

III. VORTEX LIQUID PHENOMENOLOGY 

Let us now consider the vortex liquid regime. This is the regime in which Huse and 
MajumdaJ constructed their phenomenological theory for the conductivity measurements 
of Safar et al.B They begin with a force balance equation including three terms: one propor- 
tional to the current density, J (arising from Lorentz and/or Magnus forces), one proportional 
to the average vortex velocity and hence the electric field (arising from drag and Magnus 
forces), and one proportional to the second spatial derivative of the vortex velocity and 
therefore the second spatial derivative of the electric field (arising from viscous forces): 

J^(r) = a^MKir) - S^^fs.d^dfsE^iv). (3.1) 

Although this equation was initially motivated by considering only the electric field due 
to vortex motion, it is more generally just the long-wavelength expansion of the nonlocal 
Ohm's law (1.1). The nonlocal conductivity to order k"^ is then 

cr/,^(k) = cr^^(O) + S^ap,yKkf3- (3.2) 

Note that in this paper we call the coefficient of the last term S, not rj as in Huse and 
Majumdar,El in order to avoid confusion with either the Bardeen-Stephen drag coefficient or 
the vortex liquid viscosity tensor. 

How is S related to the hydrodynamic viscosity tensor of the vortex liquid? This can 
be answered for the components of S that couple to electric fields in the xy plane using the 
more detailed force balance equation, similar to that discussed above for the lattice, given 
in Marchetti and Nelson.cl Again, B = Bz. Neglecting compressibility as well as vortex 
segments which are not parallel to the z-axis and any hexatic bond-orientational order, the 
only difference between the dc force balance equation for the liquid and the lattice is that 
the elastic force is replaced by a viscous force: 



B{J X z)i + 72n(z X v)i - -fiuvi + -r]ijki{djdkVi + d^diVu) = 0. (3.3) 

A note on notation: Bardeen and Stephen refer to the coefficient of the drag term as 77. 
Here and in Marchetti and Nelsonp the drag coefficients are 7's, and r] is the viscosity, 
which enters as the coefficient of the V^ term. Note also that this 1] involves interactions 
between vortices, as well as the connectivity and entanglement of vortex lines. If we neglect 
any dissipation that is not associated with vortex motion, following the same steps as for 
the lattice, we obtain 

C^xxO^) = —^[71 + (f^oVycxypkakp], (T^y{\<i) = -—f^[l2 " ^'OflyaxfiKkf^], CtC. (3.4) 

00-D 00-D 

Expressions for the resulting hydrodynamic contributions to S can be read off from these 
equations. 

However, the approximation of neglecting vortex segments not running parallel to the 
z-axis is inappropriate in a vortex liquid whose uniform {k = 0) resistivity parallel to the 
z-axis is nonzero. This resistivity is due to the motion of precisely those vortex segments 
that do not run parallel to the 2;-axis. Thus more generally we should consider a vortex 
liquid containing vortices running in any direction. In a large enough magnetic field and at 
low enough temperatures, only the field-induced vortices are present, and they are all nearly 
parallel to each other and to the 2;-axis. When a uniform current is applied perpendicular 
to the vortices they all move in the same direction, as in the case of the vortex lattice. 
However in a nonuniform [k = 0) current vortices in neighboring regions experience different 
forces. Unlike in the vortex lattice the resulting shears are not simply balanced by elastic 
forces, instead the vortices do continue to move past each other. However, their motion 
is slowed (relative to the case of a uniform current) by the vortex-liquid viscosity arising 
from connectivity, entanglement and other vortex- vortex interactions. Slower vortex motion 
means reduced flux-flow resistivity. Thus the fiux-fiow resisitivity decreases with increasing 
k and the conductivity increases. This corresponds to positive viscosities in the above 
hydrodynamic model and gives, for example, S^xxx > 0. 



At higher temperatures and in low or zero magnetic field, there will be thermally excited 
vortex line segments present with all orientations that outnumber the field-induced vortices. 
Parallel vortices interact repulsively, while antiparallel vortices attract each other. Thus the 
near neighbors to a given vortex segment are more likely to be antiparallel. Such nearby 
antiparallel segments can be joined into a vortex loop, or, in a film can form a bound 
vortex-antivortex pair. When a current is applied, those segments perpendicular to the 
current will feel the resulting Magnus/Lorentz force. If the current is in the x direction, 
a vortex segment parallel to positive z will be pushed in the negative y direction while a 
segment parallel to negative z will be pushed in the positive y direction. Both motions 
induce electric fields of the same sign along the x direction and thus contribute to the fiux- 
fiow resistivity. However, the relative motion of such antiparallel vortices is impeded by 
their being connected or entangled, as well as by the attractive force between antiparallel 
vortices. In a uniform {k = 0) current the antiparallel segments feel equal and opposite 
forces from the current, so this is when their motion is most impeded by connectivity, 
entanglement or other interactions. For a nonuniform {k > 0) current, the forces do not 
cancel and part of the motion generated is "center-of-mass" instead of relative, so is not 
impeded by the interactions. Thus we expect more vortex motion for k > 0, in this case 
where it is the relative motion of antiparallel vortices that dominates the resistivity. More 
vortex motion means a larger electric field, greater resistivity and lower conductivity. Thus 
here the conductivity for small k is maximal for uniform {k = 0) current and decreases with 
increasing k. This corresponds to negative viscosities in the above hydrodynamic model and 
gives for example S^xxx < 0. Thus we expect that the sign of the nonlocal effect (the S^ki^s) 
will change as one varies the field and/or temperature in the vortex liquid regime. This sign 
change has been confirmed for S^xxx in preliminary Monte Carlo simulations of a simple 
two-dimensional model superconductor.EJ 

The above phenomenological arguments deal with the long-wavelength (small k) behav- 
ior. At short wavelengths (large k) we expect the qualitative behavior is not sensitive to 
the various phase transitions and distinctions that affect the long-wavelength behavior. In 
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all the regimes where we can obtain the large-fc behavior, namely, the Meissner and vortex 
lattice phases in TDGL without fluctuations and the treatment below of the fluctuation con- 
ductivity in TDGL in the zero-field normal state, we find the conductivity decreases with 
increasing k in the Isnge-k regime. Thus we expect this remains true throughout the vortex 
liquid as well. This then suggests that when the /c-dependence at small k changes sign, 
the conductivity is changing from monotonic in k (in the higher-temperature, lower-field 
regime), to non-monotonic in k (in the lower-temperature, higher-field regime). 

IV. THE TIME-DEPENDENT GINZBURG-LANDAU EQUATIONS AND 

LINEAR RESPONSE 

In order to calculate the nonlocal conductivity of a superconductor we need to specify the 
dynamical equations of motion for the superconducting order parameter ip. We will adopt 
the simplest such description, the time-dependent Ginzburg-Landau (TDGL) equation: 

T-\dt + ij(P)ij = |;^(V - ^j^fi^ -ai^- m'i^ + C, (4.1) 

where is the scalar potential, m is the effective mass of a Cooper pair, e* = 2e is the charge 
of a Cooper pair, a{T) = aQ{T/Tc — 1) with T the temperature and T^ the zero-field critical 
temperature in the absence of the noise term, and F is the order parameter relaxation rate 
(taken to be real). Note that we are treating an isotropic superconductor. The stochastic 
noise term C(^)^) is chosen to have Gaussian white noise correlations, with the two-point 
correlation function given by 

(r(x,t)C(x',t')) = 2T-'kBT6^'\^ - x')<5(t - t'), (4.2) 

with the coefficient being determined by the fluctuation-dissipation theorem.li^ We will work 
in the limit of large Ginzburg-Landau parameter, k = A/^, where we can neglect the fluctu- 
ations in the magnetic field. Thus the vector potential. A, is static and is simply that due 
to a uniform magnetic field, and H = B. 
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As we are interested in the linear response of the system to an apphed electric field, 
we can calculate the conductivity matrix by using the Kubo formula, which expresses the 
conductivity as the Fourier transform of the current-current correlation function: 

Since we are ignoring fiuctuations in the vector potential, the current which appears in the 
Kubo formula is the supercurrent, 

he* (p*V 

Js = ^— (V^*VV^ - ^VV^*) - ^IV'I'A, (4.4) 

and the conductivity is due to superconducting fluctuations; the total conductivity is ob- 
tained by adding this contribution to the normal state conductivity. The validity of the 
Kubo formula in the context of the TDGL equations (with real F) can be demonstrated.E3 
However, when F is complex, the usual form of the Kubo formula may need to be modifled.li^ 
If we make the Gaussian approximation and neglect the cubic term in the TDGL equa- 
tion, Eq. ( [4.1| ), then the current-current correlation function factors into a product of order 
parameter correlation functions, with the result that 

^^ ' 2kBT \2miJ J 2tx J ^ ' 

{D^i - D*3)(L)^2 - Dl^) (:7o(x2,X3,^)Co(xi,X4,u;) 13^14^2 ' (4-5) 

where D^ = d^ — ie*A^/h, and 

Co(x,x» = l^e— (*-*')(V^(x,t)V^*(x',tO) 

2k T 
= _^IniG'o(x,x», (4.6) 

u 



with Go{:s.,x.';uj) the order parameter response function, which is the solution to 

Go(x,x» = 5(^)(x-x'). (4.7) 



-iT-'^u (V - i^Af + a 

2m n 



Within the Gaussian approximation we have also obtained the conductivity directly from 



the equation of motion (|4.1| ), as a check on the Kubo formula calculation. Thus for general 
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k in zero magnetic field and to order k"^ in a magnetic field we have confirmed that the 
fluctuation contribution to the dc conductivity matrix is indeed symmetric and given by 
( [4 .31 ). Note that the Hall conductivity is zero-M 

The Gaussian approximation is valid in a region well above (as defined by the Ginzburg 
criterion) the superconducting transition, where there are only small amplitude fluctuations 
of the order parameter. When fluctuations are large, the nonlinear terms become important 
and the Gaussian approximation no longer applies. It is precisely in this strong-fluctuation 
regime where a description in terms of fluctuating vortices becomes appropriate. This is 
the vortex liquid regime, for which this calculation is inappropriate. However, the low 
temperature end of these Gaussian results may exhibit some of the properties of the high 
temperature behavior of the technically more difficult vortex liquid regime. In principle, we 
could use the Hartree approximation^ to include a partial contribution from the nonlinear 
terms. However, because the Hartree approximation only renormalizes the critical temper- 
ature and thus will not change the signs or any other qualitative properties of S^ki, in this 
paper we will only focus on the Gaussian approximation. 

V. ZERO MAGNETIC FIELD 

Before tackling the technically difficult task of calculating the nonlocal conductivity 
tensor in an applied magnetic field, we will first calculate the nonlocal conductivity in zero 
magnetic field in the Gaussian approximation. In this case the system is translationally 



invariant, and after Fourier transforming the Kubo formula, Eq. ([4.5| ), we obtain 

2 fhe*\ r duo r d'^p ^, , , ,^, , , , , . 

^-W = ^M /5^/(^P.P^C-,(p + k/2..)C„(p-k/2..). (5.1) 

where the correlation function is 

Co(k, .;) = ^ -,. (5.2) 

(cu/r)2 + (^2p/2m + a) 

For d < A the integral is ultraviolet-convergent, so we do not need a cutoff. After substituting 

the correlation function into Eq. ( |5.1|) , performing the frequency integral, and scaling the 
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momenta by the correlation length ^ = h/ J2m\a\, we obtain 

a^,(k) = a(0)F^,(kO, (5.3) 

where the scaling function F^j^y{x) is normalized so that F^^ifi) = Sf^u, and where the k = 
conductivity iS'iB 



h^T ) (47r)'^/2 



a{0) = ksT { —^— e ■ (5.4) 



The calculation of the scaling functions F^^{k^) is rather complicated, and the details 
are relegated to Appendix A. The conductivity can be decomposed into transverse and 
longitudinal components 

a,M = ^mF^ikOP^. + F^kOPiX (5.5) 

where we have introduced the transverse and longitudinal projection operators: 

P^ = S,. - ^, P' = ^. (5.6) 



A steady-state electric field is purely longitudinal (it is the gradient of the scalar poten- 
tial), so only the longitudinal part of the conductivity enters in determining a dc current 
pattern. In the Gaussian approximation, the transverse and longitudinal scaling functions 
can be obtained in closed form (see Appendix A) and are plotted in Figs. 2 and 3. Both 
functions decrease monotonically with increasing k. As shown in Appendix A, for large x, 
F'^'^(x) ~ c^' {x/2)~^^~'^\ with c/ universal constants (there are logarithmic corrections 
in two dimensions). Expanding the scaling functions to 0{k'^), we obtain 

a(0)F^'^(A;O = (t(0) + S^'^P + ..., (5.7) 

with 

S^ = -1(4 - d)a{0)e, S^ = "1^(4 - d)a{0)e. (5.8) 

The full S tensor can be written compactly as 
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S^iajBv — S 5^jy6ap + -{S —S ){5^a^yp + 5^p5^a)- (5.9) 

Finally, after Fourier transforming back to real space, we find for the current 

J(x) = a(0)E(x) + ^^V X V X E(x) + S^WV ■ E(x) + .... (5.10) 

Note that S^ (which is equal to, e.g., Sxxxx) is negative here in the Gaussian approxi- 
mation, and is argued above (Sec. Ill) to remain negative in the vortex liquid in the critical 
regime just above Tc at H = 0. The general scaling form (5.5) presumably also remains 
valid in the critical regime, but with quantitatively different scaling functions. However, 
the qualitative behavior of the conductivity — that it is maximal for k = and falls off as 
a power of k for large k^ — should be the same in the nontrivial critical regimes for H = 0. 
Thus there is no sign here that anything dramatic occurs in the nonlocal conductivity at the 
Ginzburg crossover from mean-field to nontrivial critical behavior in zero magnetic field. 

VI. NONZERO MAGNETIC FIELD 

In this section we will calculate the nonlocal conductivity in the Gaussian approximation 
in a uniform applied magnetic field H = Hz. We will work in the Landau gauge A = 
{0,Hx,0). The response function is obtained as an eigenfunction expansion in a harmonic 
oscillator basis set in the x direction — the Landau-level expansion. Using Eq. ( [4.7| ), we then 
obtain for the correlation function in a mixed representation 



oo 



ry f 111 \ r,i rrrp-l V^ Un[X Xo)Un[X Xq) 

Co{x,x;ky,k„uj) = 2kBTT }^ ^ , (6.1) 

„,=0 [^/^) +^n,k. 



where the oscillator functions are 

\l/2 



""<^>=(2;;^j ^-'''"»*»(V'«). (6.2) 

with the energies 

h'^k'^ 1 

^n,k, = -TT^ + hujc{n + -) + a. (6.3) 

Zm 2 

In the above we have introduced the magnetic length Ih = {h/e*HY^'^, the orbit center 
coordinate xq = I'uky and the cyclotron frequency ujc = e*H/m. 
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A. J, E perpendicular to H 

In this subsection we calculate the nonlocal conductivity when the current and electric 
field are in the x — y plane. Given the symmetries of the system, in this geometry at order 
k"^ there are only three independent coefficients to be calculated: Syyyy, Syxxy, and Syzzy To 
obtain these it is only necessary to calculate ayy(k). Rotational symmetry about the z-axis 

gives Oxxxx '~'yyyyi '~'xyyx '-'yxxyi ^xzzx ^yzzyi GXC. aUQ Oxxyy ^yyxx K^yyyy ^yxxy ) / ^■ 

Due to the symmetries and the absence of a Hall effect, all Sijki with unpaired indices (e.g., 
Sxyyy) vaulsh lu thls calculation. 

To calculate ayy(k), we carry out the differentiations indicated in Eq. ( [4.5|) , (noting that 
Dy = dy — ieHx/h) using the correlation function in Eq. ( |6.1|) . The frequency integral is 
then easily performed; after scaling x and Xq by the magnetic length l^, and rescaling the 
external momenta in the x — y plane by Ih so that {kx, ky) = [kxln, kyln), we obtain 

ayy{k) = kBTr''^{e*f ( ) -- V Imn{kx,ky) 



oo 



, dPz_ 1 , ^. 

J-oo 27r SmiPz + kj2)en{pz - kz/2)[em{Pz + kz/2) + en{pz - kj2)\ ' 

where 

/oo __ /"OO _ 

d{xi - X2)e'''"(^i"^2-' / dxQ (xo - Xi){xq - X2)m„(x2 - Xq - ky/2) 
'OO J —oo 

XUn{xi - Xo- ky/2)Um{xi - Xq + ky/2)Um{x2 - Xq + ky/2). (6.5) 

After shifting the integration variables, this integral can be transformed into a single integral: 

^mn{kxj ky) = \A^^[kx, ky)\ , (o.oj 

where 

_ _ POO - _ _ 

A^lCkx, ky) = / dye-'^^yyu^y + ky/2)un{y - ky/2). (6.7) 



Some simplification also occurs in the p^-integral, if we first rescale pz by the zero- 
temperature correlation length ^(0) = {h /2ma<^^l'^ ., and introduce kz = kz^{0). Then 
we find 
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^yy(^) = /^T- E lA^^nihM'B^^Ch), (6.8) 

where we have introduced the thermal length A^ = ^Q/lGTr^fc^T,!!^ the flux quantum 0o = 
2'Kh/e*, and the scaled magnetic fleld h = H/Hc2{0), with Hc2{0) = 0o/27r,^^(O) the zero 
temperature critical fleld. The integral -B^^ is given by 

B^^lik) = Ah^ r g . . r ,.... .. r ,J. ... . .r ,... . ^^ (6-9) 



-oo 271 [{p + kj2f + /X^] [(p - kj2Y + /i„] [p2 + {kj2f + //(^+„)/2] 

where 

/i^ = ei^ + 2/im, eH = {T/T,)-l + h. (6.10) 

Within mean-fleld theory the transition to the flux lattice state occurs at e^ = 0. 

We now perform a long wavelength expansion of the conductivity. The expansion of A\^^ 
is most easily carried out by flrst noting that it can be written in a compact operator form 
(with p = \§- the momentum operator), 

A^^^{k^,ky) = {m\e~^''y^e~^''^^ye~^^''^e~^^''^\n) 

= e3^-^^(m|e-5(^-2^+^^P)^e-5(^-2^+^^P)|n), (6.11) 

where we have used [y,p] = i to obtain the last line. Expanding for small {k^, ky), we flnd 

(neglecting the multiplicative phase factor) 

- i - 1 - 

^mii^xjky) = {m\y\n) - i{m\f\n)k^ - -{m\yp + py\n)ky - -{m\f\n)kl 

— {m\yp^ + 2pyp + p^y\n)k1 
8 *^ 

--{m\?>fp + 2ypy + 2>pf\n)kjiy + 0{e). (6.12) 

o 

The oscillator matrix elements are tabulated in Appendix B. After squaring, we flnd 

\A^^^{k^, ky)\^ = - [n5m,n-l + (n + I) 5„i,n+l] + [{n - l)nS„,,n-2 - ^n"^ Srr,,n^l 



\2 



+ {2n + ly 5„,„ - 3{n + 1)' 6m,n+i + in + l){n + 2) 5„,„+2j ikj2) 

+ (n - l)n 5m,n-2 - n^ ^m,n-l - {u + 1^ 6m,n+l 

+ {n + l){n + 2) 6m,n+2] [ky/^Y + 0{kl k'y)- (6.13) 
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The long wavelength expansion of ^^^(fcz) is discussed in Appendix B. After substituting 
these expansions, Eqs. (|6.13|) and (piOp into the expression for the conductivity, Eq. ( |6.8D , 
and returning to conventional units, we obtain the k = conductivity 

1 2 1 






(a + 2n)i/2 (a + 2n + l)i/2 ^ (« + 2n + 2)V2 



(6.14) 

which agrees with previous results.li3nl3'EJ We have defined the scaling variable a = en/h, 
which measures the temperature distance from from the mean-field Hc2 line in units pro- 
portional to the magnetic field. Thus a runs from zero on the Hc2 line at nonzero field to 
infinity at if = in the normal state above T^. Note that the sum here (and in many of 
the results below) is a scaling function that depends only on this one parameter, a. This 
type of sum may be written in a more compact form by using the integral representation for 
the function $(2;, s, v)x^ After resummation, they are the Laplace transformations of certain 
functions, so we define 



fOO gfj 

aijiO) = a(0) (a - 1)^/' / -^ exp(-at) ^^,(0) , (6.15) 

Jo VTT 

Sijki = a(0) e{T) {a - If/^ / -^ exp(-at) S,,ki , (6.16) 

Jo \/7l 



where we have expressed the prefactor m^{0) /STcfi^T Arp in terms of a, the zero-field coherence 

length ^(T), and cr(0), the k = 0, H = conductivity at temperature T. We find 

4t-i/2 
(l + exp(— t))^ 

We also obtain: 

2 00 r / , 1 \ /o , r.\ A / , -1 \2 



mm Ih . ^ 

^yyyy ^^^^2^^,^1,1/2 ^2 ) ^^ 



(n + l)(3n + 2) 4(n+l) 



2(a + 2n)V2 (a + 2n + 1)1/2 



(n + l)(3ri + 4) (n + l)(n + 2) 



(a + 2n + 2)i/2 2(a + 2n + 4)V2 



-t -i/^(l-exp(- t)) 
;i + exp(-t))2 



'-'yyyy ~ (t _l ovr^(—+^^3 ' l,u.±o; 



s, 



yxxy 



+ 



C _ f-1/2 



12(n + lY 



E 

n=0 



3 (2n+l)2 (n+l)(lln + 10) 



4 (a + 2n)5/2 
n+l)(7n + 8) 



+ 



2(a + 2n)i/2 
(n + l)(n + 2) 



a 



2n + 1)1/2 (a + 2n + 2)i/2 2(a + 2n + 4)i/ 



(t^ - 5) exp(-4t) + 12 exp(-3t) 



(1 - exp(-2t))3 
(6t2 - 14) exp(-2t) + 12exp(-t) +t^ -5 
(l-exp(-2t))3 



(6.19) 



5, 



K(o) /e(o)^ 



s/^^y 



E 



fn + i) 



87r;i'rAT/i3/2 V 2 ; ^^ [ (a + 2n + 1)^/2 
1 1 



+4(n + 1) 



5. 



(a + 2n + 2)1/2 ^^ ^ 2n)i/2 
-12(n + 1) ((a + 2n)i/2 - 2{a + 2n + 1)^/^ + (a + 2n + 2)^/^) 
2t-3/2 [(2t + 3) exp(-2t) + (t^ - 6) exp(-t) + 3 - 2t] 



yzzy 



(6.20) 



(1 - exp(-2t))2 . 

An advantage of the integral representations is that one can determine the signs much 
more easily. We find that for t > 0, Sijki and dij are either always positive {Syy^z, 5-^^,(0), 

and CTzzy^)) OT always negative [!^yyyy, O^zzz^ ^zyyz^ '-'yzzy 1 >^zyyz~\~ 2ozzyy^ and Oyzzy~\~ 20yyzz)- 

The only exception is Syxxy which changes sign at about t=1.45. As a result, Syxxy changes 
sign at about a = 1.2 {Syxxy becomes positive for a ~ 1.2). 

The zero field limit (a —^ oo) can also be easily obtained in the integral representations. 
Simply by redefining at as x in ( |6.15|) and ( |6.16|) and then keeping only the leading order 
terms of the power series in x of the integrands, we can recover the zero field results of the 
previous section. 



B. J, E parallel to H 

In this geometry we calculate (722 (k). The calculation of cr22 closely parallels the cal- 
culation of cTyy, so we will only outline the results. First, we can express the conductivity 



19 



as 



^-(^) = .rf^ri ^ £ I^^U^.,yPi?il(^.), (6.21) 

6znn i At ^^„=o 

where 

_ _ poo _ _ _ 

^i^U^x, ky) = / dy e-'^-yum{y + ^y/2)u„(i/ - ky/2) 

= e^~^^~^y{m\e'^~^^y+~^yP^\n), (6.22) 

and 

B^k) = 16 r ^ . ^^ . . (6.23) 

Performing the long wavelength expansion with the help of the results in Appendix B, we 
have for A^^ 

I^^U^X) ky)\^ = ^rn,n + 2 [^^m,n-l " (2^ + 1)(5™,„ + {n + l)(5„_„+i] (fc^ + fcj) 

+0{klklkm). (6.24) 

Combining this with the long wavelength expansion of -B^^, we find for the conductivity 

'^^^ 327r;i2rAT/iV2 ^^ (a + 2n)3/2 ' 

2ti/2 



which again agrees with previous results,tZril3a and 

me'(O) ^ r n + 1/2 



a + 2n)3/2 
+4(n + 1) ((a + 2n)^/2 _ 2(a + 2n + 1)^/^ + (^^ ^ 2n + 2)^/^ 



~ _ t-^'' [(t^ - 2) exp(-2t) + 4 exp(-t) + t^ - 2] 

'^^^^^" (l-exp(-2t))2 ' ^^-^^^ 

3me^(0) ^ 1 



5127r^^rAT/i3/2 ^^ (a + 2n)5/2 ' 

-t3/2 

^^"^ = 4(l-exp(-2t))- ^^-^^^ 
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C. J perpendicular to H, E parallel to H 

We now need to calculate (Ty^(k). Following the same steps as in the previous two 
sections, we find for the conductivity 



O'yzi'^-) — s^ft2p/\_ /,l/2 ^ ^mnV^xiky) Aj^^{kx, ky)Bj^^{kz), (6.28) 

where 



Snh'T Ar h^^ ^^^^ 



oo 



B^^U~k ) = 4:h^ r~ — - (6 29) 

•^-oo 2n [{p + ~kj2y + ^„][(p - kj2y + finW + ikz/2y + /i(™+„)/2] 

Expanding the A^^s using the results of the previous sections, we find 

A^mniK, kyYA^^lik^, ky) = i[n6m,n-i + (2n + l)6m,n + {n+ l)6m,n+l](kx/2) 

-[-n6m,n-i + {n + l)6m,n+i]{ky/2) + 0(kl, k^, kjzy). (6.30) 

Combining this result with the long wavelength expansion of 5^^ in Appendix B, we obtain 

1 1 



me^(O) 1 - 



(n + 1) 



Jln H'n+l, 



4(n+2) / 1/2 1/2 _ 2^1/2 



t-3/2 [(t + 2) exp(-t) + t - 2] 
;i_exp(-2t))((l + exp(-t)) 



4... = ; „!r:r::::^,:i:^:^ (6.31) 



In zero field, Syyzz approaches {S^ — S'^)/2 which is positive. 

Finally, we find that a^y(A;^, ky, k^) = <7yz{-kx, ky, kz), so that Szzyy = Syyzz- 

D. General geometry with longitudinal electric field 

In a truly dc steady state the electric field must be purely longitudinal. Thus let us 
consider general wavevector k and ask what the current is in linear response to a longitudinal 
electric field. Without loss of generality, we can take k = kyj + /c^z, Ey = Eky/k and 
Ez = Ekz/k. To order k'^, the resulting current is 
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Jy — Eky[ayy{0) + SyyyyKy + [Syzzy + '^^ VyZz) k z\ / K , (6.32) 

Jz = i^Kz[(yzz\S^) + Jzzzz^z + \Ozyyz + ^^zzyy)ky\l K, (b.Oo) 

and, in the absence of a Hall effect, Jx = 0. Thus we see that in this geometry, the nonlocal 
effect here in the Gaussian approximation is always of the sign such that the conductivity 
for a longitudinal electric field is reduced as k moves away from 0. This is of the opposite 
sign from what we argue above occurs in the vortex liquid regime at lower temperatures. 

VII. DISCUSSION AND CONCLUSIONS 

In this paper we have examined the wavevector-dependent dc conductivity, o"(k), of 
a type-II superconductor in various regimes, using phenomenological arguments and the 
TDGL equation. There appear to be at least four qualitatively different regimes of behavior 
for o"(k): First, in the Meissner phase the conductivity is infinite at A; = and monotonically 
decreasing with increasing k. This behavior should also apply in the pinned vortex lattice and 
vortex glass phases, where there is no vortex motion in linear response to a uniform applied 
current. Second, in an ideal, unpinned vortex lattice phase, the conductivity is discontinuous 
at /c = 0, taking on the finite, flux-flow value at A; = due to vortex motion, but varying as 
/c"^ for small, positive /c, where there is no vortex motion. Here the conductivity is still a 
monotonically decreasing function of k for /c > 0, but is now non-monotonic when /c = is 
included. Third, for all temperatures above T^ in zero magnetic field and for sufficiently high 
temperatures in small nonzero magnetic fields, the qualitative behavior seen in the above 
Gaussian-order TDGL calculation applies. There the conductivity is finite and maximal 
for /c = and is smooth and monotonically decreasing with increasing k. As argued on a 
phenomenological level in Section HI above, this behavior should apply in the vortex liquid 
in a low-field regime where the dissipation is dominated by the spontaneous, thermally- 
excited vortices rather than the field- induced vortices. Last, phenomenological arguments 
suggest that in the vortex liquid regime at sufficiently low temperatures and high fields 
the conductivity is instead a non-monotonic function of k: At small k the conductivity 
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increases with k due to the large vortex-hquid viscosity that impedes nonuniform motion of 
the vortex hquid. However, at larger k, where vortex motion is not the dominant effect in 
determining the conductivity, the more microscopic behavior of a conductivity that decreases 
with increasing k prevails, as it does at large k in all regimes. It is this last vortex liquid 
regime that is the least accessible theoretically, because it is a strong thermal fluctuation 
regime that does not exist in a mean-field or weak-fluctuation treatment of Ginzburg-Landau 
theory. Thus the qualitative behavior described for the first three regimes can be obtained 
directly from the TDGL equations, while the theoretical support for the description of the 
last, vortex liquid regime is, at this time, purely phenomenological. 

What other experiments might be done to probe the /c-dependence of the transport 
properties? In a transport experiment one has access only to the surface of the sample. The 
voltage contacts can measure the electric field parallel to the surface and the current contacts 
can set V ■ J at the surface. For a bulk sample, this means one must rely on modeling to 
deduce what is going on inside the sample. However, in a film geometry, the entire sample 
is surface, so one could, in principle, have much more complete measurements and control. 
With modern microfabrication techniques, it seems a study that probes down to micron or 
shorter length scales should be feasible. We pose this as a future experimental challenge. 
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APPENDIX A: ZERO FIELD SCALING FUNCTION 

In this Appendix we will calculate the integrals which appear in the scaling function for 
the conductivity in zero magnetic field. Letting k = k^ be a dimensionless momentum, we 
have for the scaling function in dimension d 
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In order to simplify the integrals we use the Feynman parameterization^ll to first combine 
the first two terms in the denominator, and then once again to fold in the third term, with 
the result that 

1 

[p2 + (fc/2)2 + p . k + 1] [p2 + (k/2y - p ■ k + 1] [p2 + (k/2y + 1] 

We then substitute this result into Eq. ( |A1| ), change variables in the momentum integral to 
q = p + {2x — l)y(k/2) to eliminate the terms linear in k, and perform the (i-dimensional 
momentum integraLcll We are then left with 

F,Ak) = F^Ck)P^^ + F^Ck)^, (A3) 



k^ 



where 



F\k) = FTk) + (2 - d/2)-e I du. I dy |i ^ (,^,).(i ! ,,y)|3-.,. . (A5) 

and where we have changed variables to w = 2x — 1. The integrals on y can be performed, 
and the remaining integrals on w can be simplified by integrating by parts. We finally obtain 

FTk) - 2 /; .. [1 + (./2)^(1 - .^r-^ - (^^^«^, (A6) 

^ ^^~ {d-2) {k/2f ■ ^^^^ 

For (i = 4 we have F'^ = F^ = 1; for rf = 3 we have 

sin-H(fc/2)/7l + (^/2)21 Jl + (k/2y - 1 
^^W-^ (.72) ' (1/2)^ ■ <^'*' 
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Jl + ik/2f - 1 



while for d = 2 



Plots of these scaling functions are shown in Figs. 2 and 3. For small k these expressions 
have the following expansions for general dimension d: 

F^{k) = l-^{^-d)P + 0{¥), (A12) 

F^(^) = 1 - 1(4 - d)P + 0{¥). (A13) 

lo 

For large k and d j^ 2, we have 

F^{k) ~ cl{k/2y^^~'^\ F^{¥) ~ ci(^/2)-(^-'^), (A14) 

where the constants are functions of dimension rf, given by 

cl = 2tdw{l- w^/^-^ - -^, c^ = -^, (A15) 

•^ io ^ ' d-2' '^ d-2' ^ ' 

with cj^3 = vr — 2. For large k and d = 2, 

T.-x 21n2 ,,-, ln(fc/2) ,, , 

The large A; behavior agrees with the result of the scaling theory discussed in Sec. Ill, up 
to some logarithmic corrections in two dimensions. 

APPENDIX B: EXPANSIONS FOR INTEGRALS IN A MAGNETIC FIELD 

In this Appendix we will include some of the details of the calculation of the viscosities 
in a magnetic field. First, we simply list some of the harmonic oscillator matrix elements 
which are used to evaluate the integrals A^^^ and A^^: 
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{m\y\n) = -y= [Vn5m,n-1 + Vn+1 5m,n+l 



{m\y'^\n) 



(m|y^|n) 



^J{n - \)nbrn,n-2 + (2n + 1) bm,n + \J {ji + l)(n + 2) bm,n+2 

+ 3(n + 1)^/2 5™,„+i + y^(n + l)(n + 2)(n + 3)5™,„. 



(Bl) 
(B2) 



2^2 



1+3 



(m|p|n) = -j= \-\fnbrn,n^\ + Vn+T5m,„+i 



(m|p^|n) 



-y (n - l)n(5m,„_2 + (2n + 1) bm,n - ^ {n + l)(n + 2) 5m,„+2 



{m\yp + py\n) = i 



-\l{n- l)n5m,n-2 + ■^{n+ l){n + 2) 5m,n+2 



(B3) 
(B4) 

(B5) 
(B6) 



(m\yp^ + 2pyp + p^y\n) = 1= 

v2 



y{n - 2)(n - l)n5„,„_3 - n^^^ (5m,„_i 



-(n + 1)^/2 (5„,„+i + ^{n + l){n + 2){n + 3)5rr.,n+3 



(B7) 



4i 
(m|3^^p + 2ypy + 3pf\n) = — = 



-y(n-2)(n- l)n5m,„_3 - ^^^^ 5m,n-i 



+ (n + 1)3/=^ 5™,„+i + ,J{n + l){n + 2){n + 3)6m,n+3 



■ (B8) 



Next, we will consider the long wavelength expansion of the integrals -B£^, B^}^, -B^^. 
These can be efficiently evaluated by using the Feynman parameterization, similar to what 
was done in Appendix B. We then have for _B£^ 

3/i2 /•! , /•! , y 



B^rlCkz) = ^ / dw f dy — . 



(B9) 



Expanding for small kz, we find 



^mni^z 



{m — ny 



+ 



+ 



1/2 ' 1/2 1/2 

/im /in f^{m+n)/2 



^ ' A^(m+n)/2 ^ ^ 



1 1 



1/2 1/2 

/im /-in 



12 



(m — n)^ K^ 



/^r + /^nf-2/i 



1/2 
(m+n)/2 



(W2)^ + 0(^^ 



(BIO) 
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In a similar fashion we have 



B^lik 



1 .1 

dw / dy 
-1 Jo 



{kj2Y{l - w^y^) + h{m - n)wy + fi(m+n)/2 



3/2 



+3(fc^/2)2 / dw dy 



yj2y3 



1 Jo 



{kj2y{l - t(72l/2) + /i(^ _ n)wy + /i(m+n)/2 



5/2 



K' + z^i^ -2/i 



1/2 



4 



/l2(m-n)2^^™ ' ^" -^(m+n)/2J /,2(^_^)2 

3/2 
1 12/i(n^+„)/2 / 1/2 I l/2^ q. ^^{m+n)/2 

^ 1/2 "^ /.2r^ _^^i2'^^m "^/^" '' '^■ 



1/2 ' 1/2 
_jJ,m jJin 



J-/ ^ hi ^ ( 'YYi 'yi 1 

A*(m+n)/2 ^ '^ 

If m = n, then this becomes 



h^{m — n)" 



(^./2)2 + 0(^,^). (Bll) 



Finally, for B^^ we have 






(B12) 



^mnK'^z 



4 yijU'^L'^y 
1 / 1 1 ^ 



wy 



(fc^/2)2(l - wV) + K'l^ - 1^)'^y + At(m+n)/2 



5/2 



(^-^)nMy' /ii^/' 



^ f^l/S , ^1/2 _ 2^1/2 



{Km+o{ki). 



(B13) 



Note that this last integral is odd in (tti, n), in contrast to the first two integrals, and vanishes 
when m = n. 
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FIGURE CAPTIONS 

FIG. 1. Illustration of a two-dimensional vortex lattice in a nonuniform current. The 
magnetic field is parallel to the z-axis, normal to the film. The three bold lines running 
parallel to the t/-axis are current contacts. Equal currents are injected uniformly along 
the two outer contacts and the total current is withdrawn along the central contact. This 
produces a current density, J, that is uniform along y, but nonuniform along x, as illustrated. 
The forces on the vortex lattice due to the nonuniform current elastically strain the lattice, 
producing the shear displacements, u, shown. These displacements are also uniform along 
y but vary along x. Outside of the outer contacts, J = and u is uniform. 

FIG. 2. Zero field scaling functions at rf = 3. 

FIG. 3. Zero field scaling functions at rf = 2. 
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